arXiv:1504.04812v2 [astro-ph.EP] 9Jul2015 


Draft version July 10, 2015 

Preprint typeset using DTeX style emulateapj v. YljX^lW 


MILLIMETER-WAVE POLARIZATION OE PROTOPLANETARY DISKS DUE TO DUST SCATTERING 

Akimasa Kataoka Takayuki Muto^, Munetake Momose Takashi Tsukagoshi Misato Fukagawa^’®, 
Hiroshi ShibaP, Tomoyuki Hanawa^, Koji Murakawa®, Cornelis P Dullemond ^ 

^Institute for Theoretical Astrophysics, Heidelberg University, Albert-Ueberle-Strasse 2, 69120 Heidelberg, Germany 
^National Astronomical Observatory of Japan, Mitaka, Tokyo 181-8588, Japan 
^Department of Earth and Planetary Sciences, Tokyo Institute of Technology, Ookayama, Meguro, Tokyo 152-8550, Japan 
^Division of Liberal Arts, Kogakuin University, 1-24-2 Nishi-Shinjuku, Shinjuku, Tokyo 163-8677, Japan 
^College of Science, Ibaraki University, 2-1-1 Bunkyo, Mito, Ibaraki 310-8512, Japan 
®Graduate School of Science, Osaka University, 1-1 Machikaneyama, Toyonaka, Osaka 560-0043, Japan 
^Genter for Frontier Science, Ghiba University, 1-33 Yayoi-cho, Inage, Ghiba 263-8522, Japan and 
^Gollege of General Education, Osaka Sangyo University, 3-1-1, Nakagaito, Daito, Osaka 574-8530, Japan 

Draft version July 10, 2015 

ABSTRACT 

We present a new method to constrain the grain size in protoplanetary disks with polarization obser¬ 
vations at millimeter wavelengths. If dust grains are grown to the size comparable to the wavelengths, 
the dust grains are expected to have a large scattering opacity and thus the continuum emission is 
expected to be polarized due to self-scattering. We perform 3D radiative transfer calculations to 
estimate the polarization degree for the protoplanetary disks having radial Gaussian-like dust surface 
density distributions, which have been recently discovered. The maximum grain size is set to be 
100 /im and the observing wavelength to be 870 /im. We find that the polarization degree is as high 
as 2.5% with a subarcsec spatial resolution, which is likely to be detected with near-future ALMA 
observations. The emission is polarized due to scattering of anisotropic continuum emission. The map 
of the polarization degree shows a double peaked distribution and the polarization vectors are in the 
radial direction in the inner ring and in the azimuthal direction in the outer ring. We also find the 
wavelength dependence of the polarization degree: the polarization degree is the highest if dust grains 
have a maximum size of Umax ^ A/27r, where A is the observing wavelength. Hence, multi-wave and 
spatially resolved polarization observations toward protoplanetary disks enable ns to put a constraint 
on the grain size. The constraint on the grain size from polarization observations is independent of or 
may be even stronger than that from the opacity index. 

Subject headings: dust - polarization - protoplanetary disks 


1. INTRODUGTION 

Dust coagulation is the first step toward planet forma¬ 
tion. Observational constraints on grain size is impor¬ 
tant in understanding the first stage of planet formation. 
The grain growth in protoplanetary disks has been evi¬ 
denced by a low opacity index of protoplanetary disks at 


(e.g., Beckwith & Sargent|1991 

iWilner et al.|2000 Kita- 

mura et al.||2002t| Andrews ^ Williams||2005t IRicci et al. 

2010b|a Guilloteau et al. 12011 ), 

which suggests that dust 


N akagawajl 1 9931 jD’Alessio et al.||2001j |Draine|2U(J6|) 

However, the constraints on the gram size with opacity 

index still have uncertainties. The low opacity index at 
millimeter wavelengths can be reproduced not only by 
changing grain size, but also by cha nging chemical com¬ 
position (e .g., Pollack et al. 1994) or shape (e.g., Min 
In addition, optically thick disks may ac- 


|et al. 2005). 

count for t. 


le low spec tral index of protoplanetary disks 


(e.g., Ricci et al. 2012). 


Recently, millimeter-wave observations have shown 
that some protoplanetary disks have a strong asymmetry 
of the surface brightness distribution alt hough the gas is 
more broadly and smoothly distributed dGasassus et al. 
' 20131 der Marel et al.||20131 Fukagawa et al.||201^ 


Isella et al.|2013 Perez et ai.|2014 ). The lopsided contin¬ 


uum emission with the smooth gas distribution may be a 
signature of dust accumulation by ga s density enhance - 
ment created by such as vortex (e.g., |Lyra et al.||2009|). 
The dust accumulation accelerates the gram growth and 
thus it may be a hot spot of planet formation. Therefore, 
the dust coagulation at the hot spot should be observa- 
tionally investigated. 

In this paper, we propose a new method to constrain 
the grain size in protoplanet ary disks, which make use of 
polarization due to scattering by dust grains. If grains 
are sufficiently large, the scattering opacity is large and 
therefore there is a significant amount of scattered light 
component in dust continuum emission. In such cases, 
in tandem with the anisotropy of the continuum emis¬ 
sion, we shall show that the dust continuum emission is 
polarized. 

Polarization of early -stage protoplanetary disks has 


been detected recently dRao et al. 2014 Stephens et al 
2014 Segura-Cox et al. |2015| ). It is widely believed that 


the polarization is due to the grain alignment with mag¬ 
netic field, as analog ous to the interpretations of star 


forming reg ions (e.g., Girart et al.|2006 |2009 Hull et al. 
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2013 2014). We propose another possibility of produc 
mg polarized dust continuum emission. The proposed 
mechanism will be tested with ALMA, which has a high 
sensitivity and a high spatial resolution. 

This paper is organized as follows. In Section 
we present the theoretical background of the scattering 
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and the polarization properties of dust grains. We also 
present a toy model to illustrate that the anisotropic dis¬ 
tribution of background emission produces the polariza¬ 
tion. In Section ^ we perform radiative transfer calcu¬ 
lations of protoplanetary disk models. In Section we 
discuss detectability of the polarization due to dust scat¬ 
tering. In Section we conclude this paper. 

2. THEORETICAL UNDERSTANDING OF SCATTERING 
OF ANISOTROPIC RADIATION FIELD 

In this section, we describe the theoretical backgrounds 
for polarization due to dust scattering. To obtain the po¬ 
larization, two conditions are required. One condition is 
concerning dust grains: the dust size should be appro¬ 
priate in a way that sufficient polarization degree is pro¬ 
duced when incident radiation is scattered. The other is 
concerning the radiation field: the incident light should 
have anisotropic distr ibution. After presenting our dust 
model i n Section |2.1[ the first condition is discussed in 
Section s 12.2112.41 and the second condition is discussed in 
Section 12.51 


2.1. Dust model 

Opacity of protoplanetary disks is dominated by dust 
grains. The opacity depends on dust properties such as 
their composition, shape, and size. In this paper, we 
focus on the dependence on the grain size, which dra¬ 
matically changes the optical properties. 

Dust grains are assumed to be spherical with radius a, 
and their compositi on is fixed. The compo sition is taken 
to be the same as Kataoka et al. (2014): the mixture 
of silicate, water ice, and organics, where t he mass frac¬ 
tional abundance is set to be consistent with IPollack et al.l 


( 1994D : Csiiicate = 2.64xl0-^ Corganics = 3.53xl0-^ and 
Cice = 5.55 X 10“^. T he refractive index of astronom ical 
silicate is ta ken from iWeingartner fc Draine (2001), or- 
ganics f r om lU ollack et aL| (|1994|), and water ice from 
Warren| (|1984|). We use IVlie theory to calculate the ab¬ 
sorption and scattering mass opacities and the scattering 
matrix Zij (see Appendix for more details). Dust grains 
are assumed to have a size distribution of n{a) oc 
The minimum size is taken to be amin = 0.01 /im and we 
change the maximum size Umax in the following discus¬ 
sion. Note that the minimum grain size is small enough 
not to affect the results. 


2.2. Scattering opacity of large dust grains 

Here, we illustrate the dependence of the scattering 
efficiency on the grain size. Figure shows the absorp¬ 
tion and scattering mass opacities or dust grains which 
have the maximum sizes of 1 jam and 100 jam. The mass 
opacity is written in unit of cm^ per gram of dust grains. 
In the case of Umax = 1 /^ni, the scattering opacity is 
much less than the absorption at ^ mm wavelengths. In 
the case of a^ax = 100 by contrast, the scattering 
opacity is larger than the absorption. This implies that 
the scattering at mm wavelengths is efficient if the dust 
grains are as large as 100 fim. This is realized in proto¬ 
planetary disks because dust grains are believed to hav e 
grown to millimeter size (e.g., Beckwith fc Sargent|1991 ). 


2.3. Polarization dependence on the scattering angle 



Fig. 1. — The absorption and scattering mass opacities of dust 
grains in the cases of a-max = 1 ^^<1 o-max = 100 p,m. The size 

distribution is n{a) oc 


The polarization by dust scattering strongly depends 
on its scattering angle. We investigate the dependence of 
the polarization on the grain size, the scattering angle, 
and the wavelengths. Figure shows the degree of po¬ 
larization due to single scattering averaged with the size 
distribution (—^i2,efr/^ii,eff) for several maximum grain 
sizes in the cases of A = 3.1 mm and A = 870 pm. In 
the cases of Umax = 10 pm and 100 /rm, the polarization 
degree shows a peak at the scattering angle of 0 = 90°, 
which is the case for the Rayleigh scattering regime that 
is realized when the grain size is smaller than the wave¬ 
length. More spe cifically, the size paramet er x = 27ra/A 
is less than unity (Bohren & Huffman|1983). In the cases 
of ^^max = 1 mm and I cm, by contrast, the polariza¬ 
tion degree is 0 for almost all the scattering angle except 
at around 135°, where the polarization degree is around 
—0.2. Note that at the minimum, the sign of the polar¬ 
ization has be en changed, which is known as the polariza¬ 
tion reversal (|DanielJ[1980| [Bastien fc Menard ||1988| |Fis- 


cher et al.||19y4t |lVlurakawa||2010t |Kirchschlager Wolf 

2014p . Thus, no polarization is expected if the maximum 
gram size is larger than roughly the wavelengths, except 
for the polarization reversal. 


2.4. Detectable grain size for each wavelength 

We now discuss at what grain size the polarization due 
to scattering is significant. On one hand, we have demon¬ 
strated that the scattering by dust grains is significant 
only when the grain size exceeds roughly the wavelength 
as shown in Section |2.2[ On the other hand, the grain 
size should be less than ^ A/27r to have significant po¬ 
larization in the scattered light as shown in Section |2.3[ 
Thus, there is a grain size which contribute most to the 
polarization. 

We investigate the dependence of polarization effi¬ 
ciency on grain size especially in the cases that wave¬ 
lengths are 870 pm and 3.1 mm, which correspond to 
ALMA Band 7 and 3, respectively. Figure shows both 
the albedo uj = /^sca/(^abs + ^sca) and the polarization 
P = — ^i 2 ,eff/^ii,eff at 90°. If the grain size is small, the 
scattered light is perfectly polarized in the case of 90° 
scattering. When Umax ^ A/27r, the polarization starts 
decreasing to 0. By contrast, the albedo uj increases with 
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X=3.l mm (ALMA Band 3) X=870 |im (ALMA Band 7) 




scattering angle [degree] scattering angle [degree] 


Fig. 2.— The degree of polarization —Z 12 /Z 11 of dust grains against scattering angles. The left panel is for the case of A = 3.1 mm and 
the right panel for A = 870 p,m. The size distribution is assumed to have power law of n(a) oc a~^'^. The maximum grain sizes are 10 /rm, 
100 /rm, 1 mm, and 1 cm. 
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Fig. 3.— The polarization P at the scattering angle of 90° and the albedo uj = Ksca /(^abs + i^sca.) as a function of maximum grain size. 
The size distribution is assumed to have power law of n{a) oc a~^'^. The wavelengths are assumed to be 3.1 mm for the left panel and 
870 pm for the right panel. The arrows indicate the maximum grain size which has the most efficient polarization by 90° scattering. 


TABLE 1 

The detegtable grain size for observed 

WAVELENGTHS 


wavelengths A 

the detectable grain size Umax 

7 mm 

1 mm 

3.1 mm 

500 pm 

870 pm 

150 pm 

340 pm 

70 pm 


increasing grain size. 

The product of the polarization and the albedo, Pcj, 
gives the grain size that contributes most to the po¬ 
larization at any observing wavelength because the ac¬ 
tual observed quantities are Stokes parameters Q and 
P, which is proportional Puj in the case of single scatter¬ 
ing. In other words, Puj represents a window function for 
the grain size detectable with polarization observations. 
Figure shows Puj at the wavelengths of A = 340 /im, 
870 /im, 3.1 mm, and 7 mm as a function of ttmax- The 
grain sizes that give the maximum values of Puj for each 
observing wavelengths are summarized in Table This 
suggests that detection and non-detection of polarization 
for a wide range of sub-mm, mm, and cm wavelengths 



0.001 0.01 0.1 1 

Maximum grain size [cm] 


Fig. 4.— The polarization degree times the albedo Puj against 
the maximum grain size. This figure represents a grain size which 
contributes most to the polarized intensity. Each line corresponds 
to the wavelengths of 0.34 mm, 0.87 mm, 3.1 mm, and 7 mm. The 
band numbers correspond to the ALMA band numbers for each 
wavelength. 
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can put a strong constraint on the grain size. 


2.5. Anisotropic radiation field 

The second condition for producing the polarization 
due to scattering is the asymmetry in the distribution 
of the light source. In the case of sub-mm emission 
from protoplanetary disks, the light source is the thermal 
emission of dust particles. Therefore, its spatial distri¬ 
bution should be anisotropic to produce polarized scat¬ 
tered light. For example, if the incident radiation from 
one direction is stronger than that from the direction dif¬ 
ferent by 90°, the final scattering is partially polarized. 
Such asymmetry or anisotropy may occur in protoplane¬ 
tary disks showing ring-like or lopsided surface brightness 


(e.g., Casassns et al. T 2 M 11 Ivan der Marel et ah 


IM3 I 

4 al.llz' 


Fuk- 


agawa et al. \2m\ llsella et al.||2013l |Perez et al.pi4|) 
in this section, we demonstrate that polarization due to 
self-scattering occurs in such anisotropic radiation field 
by using a simple toy model. Note that the polariza¬ 
tion due to scattering of anisotropic radiation field has 
been discussed in the context of E-mode polarizati on of 
the cosmic microwave background radiation (e.g., Rees 
1968 Seljakjp^^ Kamionkowski et al.||1997| ). 


Hereafter, we hx the maximum grain size and wave¬ 
lengths to be Umax = 100 and A = 870 /im, which 
is one of the best combinations of the efficient polar¬ 
ization, in order to investigate possibilities of detecting 
mm-wave polarization from protoplanetary disks. Note 
that the calculated absorption and scattering opacities 
are /^abs = 0.51 cm^/g and z^sca = 1.22 cm^/g. We stress 
here that we consider the same size distribution of dust 
grains of n{a) oc in the following discussion. 

Some protoplanetary disks show that radially Gaussian 
but azimuthally smooth surface brightness (see Section 
l^for the protoplanetary disk models). To investigate a 
possible polarization in that location, we consider a tube- 
like density distribution as a toy model. The surface 
density has a Gaussian form in two directions (^, z) as 


Pd = Po exp 


exp 


( 1 ) 


around the tube at p ^ 60 AU and y ^ —60 AU, then it 
goes up to a constant value in the outside. In addition, 
polarization vectors are in the y direction at the peak of 
the surface density, p = 0, but it is in the x direction in 
the outer side, p ^ 70 AU. 

Figure illustrates the polarization mechanism 
schematically. First we consider the last scattering at 
point A (at p = 0) before the light reaches the observer 
in the 2 ; direction. If the incident photons propagate in 
the X direction before the last scattering (the thick red 
arrows), they are polarized in the y direction after the 
scattering (the thin red arrow). If the incident photons 
are in the y direction (the thick blue arrows), they are 
polarized in the x direction (the thin blue arrow). The 
former photons dominate over the latter because point A 
is at the peak surface density and therefore, the radiation 
is partially polarized in the y direction. The polariza¬ 
tion degree depends on the degree of incident radiation 
anisotropy, or, in particular, the quadrupole moment. 
Next, we consider the last scattering at point B in Fig¬ 
ure which corresponds to the periphery of the tube. 
Most of the incident photons of the last scattering come 
from the tube region {y ^ 0) and therefore are in the y 
direction. Thus, the radiation from point B is polarized 
in the x direction. 

In summary, we have shown that the polarization de¬ 
gree is significantly high when the following conditions 
are met. One is that dust grains have a maximum size 
of <^max ^ A/27r. The other is that flux coming from two 
opposite directions are stronger than that from the 90° 
different two directions (point A in Figure]^ or that the 
radiation field has a strong gradient (point B in Figure 

§• 

3. IMPLICATIONS FOR POLARIZATION OBSERVATIONS 
OF PROTOPLANETARY DISKS 

In this section, we investigate whether the polarization 
due to dust scattering is effec tive by using protoplanetary 
disk models. In Section [3Tj we discuss the possibility of 
detection of the polariz ation using an axisymmetric ring- 
shaped disk. In Section [ 3 ^ we use lopsided disk models 
to calculate the polarization. 


The adopted values are yo = zq = 27 AU and po = 
10“^^ g cm“^. The observer is in the 2 ; direction. Radia¬ 
tive transfer calculations are performed with RADMC- 
3EB We confirm the validity of RA DMG-3D by perfo rm- 
ing a benchmark test proposed by|Pinte et al ||2009 (see 
Appendix). Note that we include tne multiple scattering 
and solve the Stokes parameters for each single scatter¬ 
ing. 

Figure shows that the intensity, polarized intensity, 
and polarization degree overlaid with polarization vec¬ 


tors 


pol 

R ^ 


^ The map of the polarization degree shows that it 
is peaked around ^ = 0, which coincides with the peak of 
the surface density. The polarization degree has two dips 


^ RADMC-3D is an open code of radiative transfer calculations 
developed by Cornells Dullemond. The code is available online: 
http://www.ita.uni-heidelberg.de/~dullemond/software/radmc- 
3d/ 

^ The polarization degree is higher at x ~ —300 AU and x ~ 
300 AU, which is due to the boundary effects but it does not affect 
the discussion. Also, the polarization vectors have different lengths 
at y = 0, which is due to drawing program but this does not affect 
the discussion. 


3.1. Polarization from a ring-shaped protoplanetary disk 

We assume a protoplanetary disk where the radial pro¬ 
file of dust continuum emission is described b y Gaussian, 
which is realized in some observations (e.g., Isella et al 
2012). The grain density distribution is assumed to be 



1 , 

U - ?’d V 

1 

U ^ V( 

Pa = So exp 

2 ^ 

Ua )J 

— exp 
v27rhg 

. 2 Ug/' 


-U w (?) 


where r is the orbital radius. The mass of the central 
star is taken to be 1.9 Mq. The adopted values are 
Eo = 0.6 g cm“^, Td = 173 AU, = 27 AU, and 
hg = 19.8 AU (r/173 AU)^’^, which corresponds to the 
isothermal disk of 36 K. The dust mass of the disk is 
5.0 X 10“^ Mq. We confirm that the results in the case 
of a power-law temperature distribution do not show 
any significant difference from the constant temperature 
adopted here (see Appendix [B] for more details). 

These choices of parameters are m otivated by recent re¬ 
sult s of the modeling of HD 142527 (|Muto et al.|submit- 
ted), although we use different dust models. In addition. 
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Fig. 5.— The intensity, polarized intensity, and polarization degree with polarization vectors in the case of the tube model. 


y 

A: at the peak 


y=0 



B: outside 





zO- 

Fig. 6.— A schematic picture of the mechanism of the polariza¬ 
tion due to scattering of anisotropic radiation field. 

the dust density is assumed to be zero if i? < 70 AU or if 
R > 300 AU. We assume that the target is at 140 pc and 
thus 1 arcsec = 140 AU. Note that the optical depth at 
the peak is r = Sq x K^hs = 0-b g cm“^ x 0.51 cm^ g“^ = 
0.31. Thus, this object is optically thin. 

Figure shows the intensity, the polarized intensity, 
and the polarization degree overlaid with polarization 
vectors. The polarization degree has a double-ring struc¬ 
ture. The polarization vectors are orientated to totally 
opposite directions in the two rings. The vectors in the 
outer polarization ring are in the azimuthal direction. 
This is because the background thermal emission has a 
strong radial gradient at the location of the outer po¬ 
larized ring. This corresponds to the point B in Figure 

t By contrast, the vectors are in the radial direction in 
e inner polarization ring. This is due to the net flux 
from the azimuthal direction is larger than the net flux 
in the radial direction. This corresponds to the point A 
in Figure This double-ring pattern is a unique feature 
of the polarization due to dust scattering, and thus this 
will be a clue to distinguish the polarization mechanism. 

3.2. Polarization from lopsided protoplanetary disks 

We now calculate the expected polarization from a 
lopsided disk. To mimic the lopsided disk structure 


observed w ith the ALMA observations (e.g., Fukagawa 
et al.|2013| ), we further add an az imnthally Gaussian dis¬ 
tribution (e.g., Perez et al.||2Q14| ) as follows: 


Pd = 


:exp 


{0 - Od) 

0d 


2i 




( 3 ) 


\/^hc 


exp 


-- f—Y 

- 

. 2 Ug/. 

exp 


r-rd 

Wd 


Since broad range of azimuthal contrast of dust con¬ 
tinuum emission has been reported so far (1.5 for 

24 for HD 142527; [FuF 


SA02Q6462; 


agawa et al. 


Perez et al.[ 

20131 130 foiMRg 


^ 48; |van der Mar el et al. 

2013p , we consider two cases: model A for low azimuthal 
contrast and model B for high azimuthal contrast. In 
the model A, we use the model where the disk is entirely 
optically thin, while in the model B, the disk is partially 
optically thick, and therefore, we can also investigate how 
the difference of optical depth affects the polarization. 

Specific parameters are as follows. In model A, the 
ratio of the surface densities between the top and the 
bottom is 4. The adopted values are Umax = 0.6 g cm“^, 
^min = 0.2 g cm“^, Td = 173 AU, Wd = 27 AU, hg = 
19.8 AU (r/173 AU)^’^, Od = 0, and 0d = 7r/3. The dust 
mass of the disk is 3.7 x 10“^ Mq. The radial peak of 
the dust surface density is located at r = r^, and it is 
maximized at P = and minimized at d = + tt. We 

call the location of the former peak top peak and the 
latter bottom peak. The optical depth at the top peak is 
T = (Smax + Smin) X Kabs = 0.8 g Cin“2 x 0.51 CHl^ g”! = 
0.41 , and the optical depth at the bottom peak is r = 
0.21 g cm“^ X 0.51 cm^ g“^ = 0.11. Thus, model A is 
showing totally optically thin case. 

In model B, the ratio of the surface densities between 
the top and bottom is taken to be roughly 70. This kind 
of the str ong enhancement is re alized in HD 142527, for 
example (Muto et al. submitted). The adopted values are 


Smax = 3.3 g cm“^, Emin = 0.046 g cm“^, and the other 
values are the same as model A. The dust mass of the disk 
is 1.2 X 10“^ Mq. The difference between model A and 
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Fig. 7.— The same as Figure but in the case of the ring-shaped protoplanetary disk. 


B is the optical depth. In model B, the optical depth at 
the top peak is t — (^max ^min) ^ ^abs — 3.3 g cm X 
0.51 cm^ = 1.68, which is optically thick. 

Figure shows the intensity, polarized intensity, and 
polarization degree with polarization vectors of the lop¬ 
sided disk model A. The polarization degree mainly 
shows a double-peaked structure. The polarization de¬ 
gree is higher in the brighter region. The polarization 
vectors are in the azimuthal direction in the outer ring, 
which corresponds to the point B in Figure and are in 
the radial direction in the inner ring, which corresponds 
to the point A in Figure This feature is the same as 
the ring model. 

Figure is the same as Figure but for model B. The 
polarization degree maps of the two models are qualita¬ 
tively similar except for that the polarization degree is 
lowered in the region of the intensity peak in model B. 
This feature is realized due to the difference in the polar¬ 
ization degree between optically thick and thin regions. 
The polarization degree at a specific location of the disk 
is determined by the anisotropy of the incident radiation 
field before the last scattering, which corresponds to the 
region within the range of r < 1 around that location. 
If the disk is optically thick, the geometric distance cor¬ 
responding to r = 1 is small and therefore the incident 
radiation field is nearly isotropic. This is the reason why 
the polarization degree decreases at the optically thick 
region. This reduction of the polarization degree results 
in two peaks of the polarization degree in the inner ring 
of model B, which correspond to the location of r ^ 1. 
This case can be tested observationally with HD 142527 
for example. 


4. DISCUSSION 

4.1. Required spatial resolution and sensitivity to deteet 
the polarization 

We have seen that the ^ 2.5 % of polarization degree is 
expected from the dust continuum emission from proto¬ 
planetary disks in which the dust particles are grown to 
sufficiently large size and are distributed in a ring-like or 
lopsided shape. Here, we discuss the required spatial res¬ 


olution to detect the polarization due to dust scattering 
with actual observations. 

First, we consider the case that we observe a disk with¬ 
out spatially resolving the entire disk. In the case of the 
ring disk model, no polarization is expected because the 
polarization is totally axisymmetric and the polarized 
flux is cancelled out. In the case of the lopsided disks, by 
contrast, the net polarization pattern is expected even 
without spatially resolving the whole disk because the 
polarization pattern is not axisymmetric. However, the 
net polarization degree is ^ 0.1 % for model A and ^ 0.06 
% for model B, respectively. Such low polarization de¬ 
gree is unlikely to be detected with instruments to date. 
Therefore, without spatially resolving the entire disk, no 
polarization is expected. 

Then, how small should be the spatial resolution to de¬ 
tect the polarization? Figure shows the Stokes Q and 
U images of the ring disk mood and the model A of the 
lopsided disk model. In the case of the ring disk model, 
the Stokes Q and U maps show butterfly-like signatures. 
To detect the polarization, we have to resolve the two 
positive peaks and the two negative peaks. Especially in 
this case, the peaks have a length of ^ 1 arcsec in the 
azimuthal direction and ^0.3 arcsec in the radial direc¬ 
tion. Thus, the spatial resolution should be high enough 
to safely resolve a 0.3 arcsec structure. Moreover, the 
dust continuum emission should be detected with good 
signal-to-noise ratio in order to detect the Stokes Q and 
U. Since the polarization degree is at most ^ 2 %, the 
Stokes I should be detected with at least ^ 250 <j level to 
detect Stokes Q and U with 5cr level. Such requirements 
of high spatial resolution and good sensitivity can be met 
in near future ALMA observations. The same discussion 
is applicable to the case of the lopsided disk model. 


4.2. Grain size eonstraints 

The polarization due to self-scattering is an indepen¬ 
dent way to constrain the grain size fro m the opacity in¬ 
dex P (e.g., jBeckwith & Sargent || 199 1|). Although 1 
at millimeter wavelengths has widely been regarded as 
an indication of dust growth, in the case of HD 142527 
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Fig. 8.— The same as Figurej^but for model A in the case of the lopsided protoplanetary disk. The object is optically thin everywhere. 
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Fig. 9.— The same as Fig. [^but for model B. At the top peak, the optical depth exceeds unity, which is realized in HD 142527, for 
example. Note that the polarization degree is saturated at the outer ring in the right figure to emphasize the polarization degree of the 
inner ring. The outer ring is unlikely to be detected because the polarized intensity is weak. 


for example, the small (3 can also be reproduced by dust 


model of ^ 3 iim. in size with an ir regular shape (e.g., Min 


et al. 2005 Verhoeff et ah 2011). However, the irregu- 


lar shape grams are not responsible for the polarization 
because of their small scattering opacity at millimeter 
wavelengths due to its small size. Therefore, polariza¬ 
tion observations clearly determine whether the grains 
are large or small with irregularly shape. 

In addition, the dust polarization is a unique and pow¬ 
erful way of determining the grain size at optically thick 
regions. The spectral index at millimeter wavelengths 
has no information of dust grains if the emission is opti¬ 
cally thick. By contrast, although the polarization degree 
has a maximum value at r ^ 1, the polarization can be 
still detectable even at the optically thick regions (see 
Figure [^. Therefore, the polarization due to dust scat¬ 


tering is a promising indicator of grain size at shorter 
wavelengths or at the inner region of disks where the 
emission is expected to be optically thick. 

4.3. Other polarization mechanisms 

We have shown that the dust scattering can produce 
the polarization of the dust continuum emission. This 
is a totally different mechanism from commonly used in¬ 
terpretation of the polarization as coming from the dust 
alignment. Our models indicate some observational char¬ 
acteristics. First of all, our mechanism requires the asym¬ 
metry of the radiation field, which comes from, for exam¬ 
ple, ring-like or lopsided distribution of continuum emis¬ 
sion. Secondly, the polarization pattern shows multiple¬ 
ring structures. Thirdly, the polarization vector is differ¬ 
ent by 90° at the two polarization rings. These are unique 
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Fig. 10.— The Stokes Q and U maps. The upper panels are in the case of the ring-shaped disk model and the lower panels are in the 
case of the model A of the lopsided disk model. 


features of the polarization due to dust scattering, and 
may be a clue to distinguish this mechanism from dust 
alignment. For example, if the polarization of disks is 
due to the grain alignment with the magnetic field, there 
is a prediction that the polarization degree is ^ 2% and 
the polarization vectors are directed in radial direction if 
the magnetic field is do minated by toroidal components 
(Cho & Lazarian||2007). This characteristics is basically 
the same as the self-scattering mechanism. Therefore, 
it is still difficult to distinguish between the two mecha¬ 
nisms and thus further and detailed comparisons should 
be necessary. 


4.4. Effects of dust settling 

We also investigate the effects of dust settling. If 
dust grains are large enough to be decoupled from the 
gas, they settle to the midplane. However, dust grains 
are stirred up by the gas turbulence. Then, the ver¬ 
tical distribution of the dust grains are approximated 
by the Gaussian function with a dust scale height, hd, 
which is determined by the balance of settling and stir¬ 
ring. The dust scale height is related to the gas scale 


height as hd = ^/a/Sthg, where a is the viscous param 
eter, St is the Stokes number of the dust grains, and hg 


is the scale height of the gas (e.g., Youdin & Lithwick 
2007). The Stokes number is given by St ^ apint/^gas/ 


where pint is the internal density of the dust grains and 
Ugas is the gas surface density. The equation is valid 
if St ^ 1 and a < St. In this paper, the grain size 
is taken to be 100 pm, which has a Stokes number of 
St ^ 10“^ X (Sgas/10 g cm“^)“^. If we assume the gas 
surface density of Ugas = 10 g cm“^, dust settling is ef¬ 
fective if a < 10“^. Therefore, the dust settling is im¬ 
portant only if the turbulence of the gas disk is as weak 
as (a ^ 10“^. 


Although the settling is not expected to occur in this 
case, it is important to investigate the effects of dust set- 
tlin g. We did the same simulations discussed in Section 
|3.1| except for the scale height, which we artificially re¬ 
duce with a factor of 3. Figure shows the results of 
the geometrically thin disk. The polarization degree is 
1.2 % in the inner ring of the polarization map although 
it is 2.2 % in the case of no settling (see Figure]^. Thus, 




















































Polarization of protoplanetary disks 


9 


the dust settling reduces the polarization degree. 

The physical reason is as follows. Although the optical 
depth in the line of sight does not change if the scale 
height changes, the optical depth perpendicular to the 
line of sight increases because the dust density increases. 
Thus, the incoming flux to be scattered by dust grains be¬ 
comes more isotropic with decreasing scale height. This 
is the reason of the less polarization degree in the case of 
the thinner disk. 

5. SUMMARY AND OUTLOOK 

We have presented a new mechanism of millimeter- 
wave polarization of protoplanet ary disks, which is due 
to dust scattering of anisotropic radiation field, result¬ 
ing from a ring-like or lopsided dust continuum emission. 
Our main findings are as follows: 


• The polarization degree due to dust scattering is 
the highest at Umax ^ A/27r, where Umax is the max¬ 
imum grain size with a size distribution of n{a) oc 
and A is the observing wavelengths. This is 
because dust grains should be sufficiently large to 
have a high albedo but should be sufficiently small 
to be in the Rayleigh scattering regime in order to 
make the emission polarized. The grain sizes that 
contribute most to the polarization are shown in 
Figure and Table The strong dependence of 
the dust size on polarization implies that we can 
use the polarization signature as an indicator of 
dust size in protoplanetary disk. 


• Even without a specific light source, the contin¬ 
uum emission is polarized due to self-scattering if 
the dust emission exhibits quadrupole anisotropy 
around each location of the disk (see Figure]^. To 
demonstrate the self-scattering, we used a tul^-like 
density distribution and performed radiative trans¬ 
fer calculations. As a result, the polarization vec¬ 
tors are perpendicular to the tube at the peak due 
to the strong flux from the tube. By contrast, the 
polarization vectors are parallel to the tube at the 
periphery of the tube because the net flux is in the 
direction of the gradient of the surface brightness. 


• In the case of a ring disk which has a radial surface 
density distribution expressed by a Gaussian func¬ 
tion, the polarization is expected due to the dust 
scattering. The polarization map shows a double¬ 
ring structure. The polarization vectors are in the 
azimuthal direction in the outer disk and in the ra¬ 
dial direction in the inner disk. The polarization 
degree has a peak value just inside the peak of the 
Stokes I emission. The polarization degree is up to 
^ 2.5% in the specific case. 


• We also performed radiative transfer calculations 
on lopsided protoplanetary disks. We use two mod¬ 
els: totally optically thin model (model A) and 
the partially optically thick model (model B). The 
model B is motivated by recent millimeter-wave 
observations on HD 142527. They also show the 
double-ring structure of the polarization. In the 
case of the model B, the emission is polarized even 
at the optically thick region. In addition, the peak 
polarization is 2.5%, which is likely to be detected 
with ALMA observations. The polarization degree 
decreases if the dust is settled because the optical 
depth increases with decreasing dust scale height, 
which makes the incoming flux more isotropic. 


We have shown the results of calculations mainly in 
the case of Umax = 100 jiim and A = 870 /im as an il¬ 
lustrative example. This set of parameters is one of the 
most efficient case for producing the polarized light. The 
maximum grain size that is responsible for producing po¬ 
larization is given by Umax ^ A/27r (see also Figure]^ 
so we expect that we can apply our results to observa¬ 
tions at different wavelengths, but further study should 
be necessary for observations at different wavelengths. 
In addition, we have assumed that the disk is face-on for 
simplicity. Since the polarization degree depends on the 
scattering angle, which is determined by the disk inclina¬ 
tion and the position angle of the major axis of the disk, 
further parameter study is necessary to apply our results 
to actual observations. 

The polarization due to scattering ha s a potential to 
determ ine the porosity of dust grains. [Kataoka et al.| 
(2014) have shown that to constrain the porosity of dust 
aggregates, we need both information of the absorption 
and scattering opacities. The continuum and polarized 
intensity have the information of both opacities. There¬ 
fore, the proposed mechanism should be further tested 
with including the porosity to constrain the porosity of 
dust aggregates in protoplanetary disks. This is a good 
observational test on the planet esimal formation scen ario 
via fluffy dust aggregates (e.g ., Kataoka et al.||2013). 

As discussed in Section [T^ the polarization is a pow¬ 
erful technique to constrain the grain size compared with 
the opacity index. This work suggests that spatially re¬ 
solved polarization observations toward protoplanetary 
disks opens a new window on grain size constraints in 
protoplanet ary disks. 
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(27 • 213) This work is supported by MEXT KAKENHI 
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Fig. 11.— To check the effects of grain settling, we artificially reduce the scale height of the disk a,s H = Hg/3. This results in less 
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APPENDIX 


SCATTERING MATRIX 

The scattering properties of dust grains are described by the scattering matrix Zij . We briefly review the scattering 
matrix to understand which matrix element determines the polarization due to scattering. Let us consider the case 
when an incident light, which has Stokes parameters {li^Qi^Ui^Vi)^ is scattered by a single dust grain. We assume 
that the d ust grains are spherical grains, and then the independent elements are reduced to 4 elements (see Bohren & 
e.g., |1983|). Choosing certain observational axes, the Stokes parameters of the scattered light (Jg, Ls) 

are given by 

/ Z\\ Z\2 0 0 



m, 


gram 


(P 


^12 ^11 0 0 

0 0 Z 33 Z 34 

\ 0 0 —^34 ^33 


ik \ 

Qi 

\y^J 


(Al) 


where d is the length from the observer, and rUgrain the grain mass. We use the scattering in atrix of which 
follows the notations of the RADMC-3D. The traditional notation Sij ( Bohren fc Hnffman||1983 ) is related as Zij = 
5'ij/(/c^mgrain): whcrc k is the wave number. Note that Zij is a function of scattering angle and dust size. 

We derive the effective scattering matrix with a power-law size distribution of dust grains. To address the dependence 
of the scattering matrix on the grain size a and the scattering angle 0, we write Zij = Zij{a^ 0). 

The total scattering opacity /^sca(<^) is given by the integration of Zii(a,6>) over the solid angle as 


^sca(^) — ^ ((2, 


(A2) 


by 


If we consider a size distribution of dust grains with n{a) (x a the total scattering mass opacity /^sca,total is given 


^sca,tot — 


Pd Jamin 


4«) 


a (2a 


'scaV^/''^grain 


{a)da 


Substituting Eq.(A2), we obtain 


where 


/ Q-max pCLmax 

i^sca{cL)ci^~^da/ / a^~^da. 

- min Q-min 

rCLma-x. / C \ 

:a,tot= / i (b Zii{a, 0 )dQj a^~‘^da/ / a^~‘‘da 

/ / PCbmax PQ-max \ 

(/ Zii{a,e)a^-‘‘da/ a^-‘>da] d^ 

\’J CLjnin Q,min / 

pO'max pO'niax 

Zu,es{ 0 )= / Zu{a, 0 )a^~‘^da/ / a^-^^da. 

'J Q-min Q-min 


(A3) 


(A4) 


(A5) 


We use Mie theory to calculate the effective phase function Zii^eff(^) in this paper. The other scattering matrix 
elements Z^j^eff(^) is derived as the same way. 

If the inci(lent li ght is unpolarized, that is Qi = Ui = Vi = 0, the Stokes parameters of the scattered light can be 
obtained from Eq.( |Al[ ) as 

m_■ „ 

(A6) 


T _ nigrain ^ j 
2.S - ,r, 




Qs = 


d^ 


-Zi2li, 


Us = 0, and Vg = 0. Thus, the degree of polarization is given by 


P = 


VQ^ + u! 


Zi2 


^11 


(A7) 


(A8) 


Therefore, the ratio of Z12 to Zn represents the degree of polarization. 
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TEMPERATURE STRUCTURE 


To check the validity of the assumption of the constant temperature throughout the disk, we change the temperature 

shows the results in the case of the power-law temperature 
Thus, the assumption of the constant 


12 


distribution as T = 36 K x (i?/173AU)“^-^. Figure 

distribution. We do not see significant change in tKe" polarization degree, 
temperature does not significantly affect the results in the case of ring-like disks. 


I [mJy/arcsec^] PI [mJy/arcsec^] P[%] 



- 2-1012 - 2-1012 - 2-1012 

[arcsec] [arcsec] [arcsec] 


Eig. 12.— The same as Eig. |^but in the case of the power-law temperature distribution as T = 36 K x (R/173AU) 
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BENCHMARK TEST OE RADMC-3D 


To test the reliability of the radiative transfer code RADMC-3D, we p erform benchmark tests especially focusing 
on the polarization. We use a sophistic ated numerical mo del propose d by [Piute et al.| 2009 and c ompared the r esults 
with other codes, which are MCFOST (Pinte et ah 2006), MCMAX (Min et ah 2009), TORUS (Harries 2000), and 
Pinball ( Watson fc Henne;^|200l| P 


The problem and numerieal settings 

The problem is set up as follows. The dust particles of the disk has a surface density profile of = Xo(r/ro)“^‘^ with 
ro = 100 AU. Xo is determined by the disk mass and is discussed later. The vertical profile is pd = po exp(—^^/2h(r)^). 
The scale height is parametrized as h{r) = ho(r/ro)^‘^^^ with Hq = 10 AU. The inner and outer edges of the disk are 
Tin = 0.1 AU and rout = 400 AU. The dust opacity is 1 pm silicate particles with a particle density of 3.5 g cm“^, 
which is different from the opacity model of the main discussion of this paper. The refractive index is taken from 
[Weingartner fc: Draine||2001 

The only parameter is the disk mass Mdisk, and we use the disk mass of M^isk = 3 x 10“^ Mq to compare the 
results of the polarization with the other radiative transfer codes. The surface density at 100 AU, Xq, is calculated 
from Mdisk = 47rXo(100 AU)^-^[rQ(j^ — To avoid photons from being trapped at optically thick inner region of 

the disk, we truncate the disk with a radius of 3 AU. We use spherical coordinate (r, 6>,0). The radial coordinate r 
is logarithmically divided into 512 grids between 0.1 AU and 1024 AU. The polar angle 0 is linearly divided into 512 
grids in0 + e<^<7 r — e, where e = 0.01, which is set to avoid a numerical error exactly on the polar axis. The 
azimuthal angle 0 is linearly divided into 512 grids in 0 < 0 < 27r. First, we solve the temperature structure with 
Monte Carlo radiative transfer simulations and then simulate the intensity and polarization with 10^ photons. The 
calculation takes 60 hours with a CPU with 3.20 GHz. The viewing grid is set to be 251 pixel x 251 pixel for the 
square of 900 AU x 900 AU; the pixel scale is 3.6 AU. 


Results 

Figure shows that intensity maps with the same set up described above calculated with RADMC-3D, MC¬ 
FOST, MCMAX, and TORU^ The inclination angle is taken to be 69.5°. The intensity is written in the unit of 
W m“^ pixeU^, where the distance is assumed to be 140 pc. The intensity of the different codes are quantitatively 
match well as shown in FigureTo quantitatively investigate the difference, we also plot the intensities and their rel¬ 
ative differences from the average in Figurej^along the three lines on the intensity maps, which is (x, y) = (x, 180 AU), 
(x, —180 AU), and (—108 AU,^). Note that 108 AU corresponds to 30 pixels and 180 AU corresponds to 50 pixels. 
We take the average without the results of TORUS due to the low signal-to-noise ratio. The large error at large or 
small X or ^ in Figures (d), (e), and (f) is because there is no data points or a number of photons is too low. In 
Figures [m] (d) and (e), the error is statistical is within 15 %. In Figure|l4| (f), the error is larger beyond around 300 
AU, which is due to a low number of photons in this inclination set up. We should be careful to this point to discuss 
the difference of the polarization degree in the next discussion. 

Figureshows that the maps of the polarization degree with the same set up. We do not plot the results of TORUS 
due to thelow signal-to-noise ratio. Instead, we plot the results of Pinball. The maps of the polarization degree have 
a pattern due to the resonance of the phase function. The results are quite similar, but it has some differences in the 
polarization degree. To quantitatively investigate the difference, we also plot the polarization degrees and their relative 
differences from the average in Figure along the three lines that is the same as the previous discussion. Note that 
we take the average without the results of Pinball. Figures 16 (a), (b), and (c) show the polarization degree along the 
three lines, and they have a good agreement in their shape that corresponds to the resonance. The relative differences 
in Figures (d), (e), (f) are low where the polarization degree is high, which is for example at the regions around 


(250 AU < X < 350 AU) in Figure 16 (d) although the difference is as large as 1.5 where the polarization degree is 


low. This is due to the relative difference is enhanced in the polarization degree when the degree is low because the 
average itself is low. In addition, the large relative difference around y < —200 AU and y > 200 AU on the line of 
{x^y) = (—108 AU, ^) is due to a small number of photons (see Figure pH] (f) and [l^ (f)). The polarized intensity is 
mainly coming from the single scattering on the surface of the disk, and thus the photon noise is larger for increasing 
of the distance from the central star. In addition, the disk is optically thick, and thus the emission of the central star 
hardly reaches the midplane, which also makes a large error at the midplane. 

Figures pT| \T^ and show the results of the same benchmark test shown in Figures [TS} [TKl and pT| in the 
case of the inclination angle of i = 87.1°. Note that the three lines for the cuts are taken at = (x, 108 AU), 

(x, —108 AU), and ( — 180 AU,^). The results of this case also show the similar results of the previous discussion. 

The relative standard deviation in intensity and polarization is summarized in Tables and The relative standard 
deviation is calculated by dividing the standard deviation by the average value. To avoid the mw signal-to-noise ratio 
regions, the relative standard deviation is calculated in the range of —360 AU < x < 360 AU for the cases of the 
horizontal lines and —180 AU < y < 180 AU for the vertical lines. The relative difference is up to around 10 % in 
intensity. On the other hand, the relative difference is up to 30 % in polarization degree. 


3 The results of the other codes are obtained from the webpate: 4 do not plot the results of Pinball because the intensity data 

http://ipag.osug.fr/ pmtec/beuchmark/ provided. 
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TABLE 2 

The relative standard deviation in intensity 


codes 

{x, 180 AU) 

i = 69.5° 

(x,-180 AU) 

(-108 AU,?/) 

(x, 108 AU) 

i = 87.1° 

(x,-108 AU) 

(-180 AU,?/) 

RADMC-3D 

3.7 % 

9.3 % 

10.1 % 

3.1 % 

4.2 % 

10.9 % 

MCMAX 

4.0 % 

3.4 % 

5.2 % 

5.9 % 

1.9 % 

9.4 % 

MCFOST 

4.4 % 

6.8 % 

10.0 % 

4.0 % 

4.1 % 

6.1 % 


TABLE 3 

The relative standard deviation in polarization 

codes 

(x, 180 AU) 

i = 69.5° 

(x,-180 AU) 

(-108 AU,?/) 

(x, 108 AU) 

i = 87.1° 

(x,-108 AU) 

(-180 AU,y) 

RADMC-3D 

11.3 % 

26.0 % 

23.9 % 

13.0 % 

11.6 % 

22.8 % 

MCMAX 

10.2 % 

12.5 % 

17.6 % 

10.6 % 

10.3 % 

16.9 % 

MCFOST 

12.0 % 

18.1 % 

27.7 % 

13.2 % 

18.1 % 

29.3 % 


From this results, we conclude the benchmark test as follows. The calculations of RADMC-3D provides the intensity 
which is quite similar to the other results. The error is statistical and the relative standard deviation is within ^ 10% 
if the number of photons is high enough. In terms of the polarization degree, the results of the different codes are 
qualitatively similar and the relative standard deviation is up to ^ 30%. If the polarization degree is low, however, 
the difference is large, which is due to the combinations of the low degree of polarization, a low number of photons 
at distant region from the star and, the optically thick disk. Therefore, it requires more photons to simulate if we 
investigate the region where the low degree of polarization. Note that the we confirm that we perform the simulations 
with enough number of photons in this paper. 
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Fig. 13.— Intensity maps of other models: (a) RADMC-3D, (b) MCFOST, (c) MCMAX, and (d) TORUS. The inclination angle is set 
to be 69.5°. 


(a) (x,y) = (x,180 AU) 


(b) (x,y) = (x,-180 AU) 


(c) (x,y) = (-108AU,y) 





(d) (x,y) = (x,180 AU) 


(e) (x,y) = (x,-180 AU) 


(f) (x,y) = (-108 AU,y) 





Fig. 14. — The intensity along the three cuts of (x, y) = {x, 180 AU), {x, —180 AU), and (—108 AU, y) in the case oi i = 70°. The lower 
panels show the difference against the average. Due to the low signal-to-noise ratio, we take the average without the results of TORUS. 
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Fig. 15. — Polarization maps of other models: RADMC-3D, MCFOST, MCMAX, and Pinball. The inclination is 69.5°. 


(a) (x,y) = (x, 180 AU) (b) (x,y) = (x,-180 AU) (c) (x,y) = (-108 AU,y) 





(d) (x,y) = (x,180 AU) 


(e) (x,y) = (x,-180AU) 


(f) (x,y) = (-108 AU,y) 





Fig. 16.— The polarization degree along the three cuts of (x, y) = (x, 180 AU), {x^ —180 AU), and (—108 AU, y) in the case oi i = 69.5°. 
The lower panels show the difference against the average. Due to the low signal-to-noise ratio, we take the average without the results of 
Pinball. 
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(a) RADMC-3D : I (i=87.1°) 
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Fig. 17.— Intensity maps of other models: (a) RADMC-3D, (b) MCFOST, (c) MCMAX, and (d) TORUS. The inclination angle is set 
to be 87.1°. 


(a) (x,y) = (x,108 AU) 


(b) (x,y) = (x,-108 AU) 


(c) (x,y) = (-180AU,y) 





(d) (x,y) = (x,108 AU) 


(e) (x,y) = (x,-108 AU) 


(f) (x,y) = (-180AU,y) 





Fig. 18. — The intensity along the three cuts of (x, y) = {x^ 108 AU), (a:, —108 AU), and (—180 AU, y) in the case of z = 87°. The lower 
panels show the difference against the average. Due to the low signal-to-noise ratio, we take the average without the results of TORUS. 
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Fig. 19. — Polarization maps of other models: RADMC-3D, MCFOST, MCMAX, and Pinball. The inclination is 87.1°. 


(a) (x,y) = (x,108 AU) (b) (x,y) = (x,-108 AU) (c) (x,y) = (-180 AU,y) 





(d) (x,y) = (x,108 AU) 


(e) (x,y) = (x,-108 AU) 


(f) (x,y) = (-180AU,y) 





Fig. 20.— The polarization degree along the three cuts of (x, y) = (x, 108 AU), {x^ —108 AU), and (—180 AU, y) in the case oi i = 87.1°. 
The lower panels show the difference against the average. Due to the low signal-to-noise ratio, we take the average without the results of 
Pinball. 











































